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The three-dimensional Anderson model with a rectangular distribution of site disorder displays two 
distinct localization-delocalization transitions, against varying disorder intensity, for a relatively 
narrow range of Fermi energies. Such transitions are studied through the calculation of localization 
lengths of quasi- one-dimensional systems by transfer-matrix methods, and their analysis by finite- 
size scaling techniques. For the transition at higher disorder we find the localization-length exponent 
v = 1.60(5) and the limiting scaled localization-length amplitude Ao = 0.57(1), strongly suggesting 
universality with the transition at the band centre, for which currently accepted values are v — 
1.57(2) and Ao = 0.576(2). For the lower (reentrant) transition, we estimate v = 1.55(15) and 
Ao = 0.55(5), still compatible with universality but much less precise, partly owing to significant 
finite-size corrections. 



I. INTRODUCTION 

According to Anderson's original scenarioS, if one 
starts from a translationally invariant model Hamilto- 
nian for noninteracting electrons and considers increasing 
quenched disorder (e.g. by assigning random self-energies 
to lattice sites), a critical threshold is reached at which 
quantum interference inhibits propagation, and the sys- 
tem undergoes a localization transition, from a diffusive 
(metallic) phase to an insulating (localized) regime. For 
the past twenty years, the framework for studies of criti- 
cal properties of the Anderson transition has mostly been 
provided by the scaling theory of localizatio: 
ical studies have been progressively refined! 
accurate estimates exist for several quantities of interest. 

In the site-disordered version which will concern us 
here, the problem is usually formulated in terms of a 
tight-binding Hamiltonian on a regular lattice, 
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where the site self-energies Si are independent, identically 
distributed random variables obeying a specified distri- 
bution, denotes nearest neighbours and the energy 
scale is set by the hopping matrix element, V = 1. The 
k> , intensity of disorder is generically represented by a width, 
j_j ■ W, of the self-energy probability distribution. 

On a simple cubic lattice, for the non-random case 
(all £j = So = 0, for convenience) the Schrodinger equa- 
tion Ti. = E \tp) yields Bloch (propagating) states for 
-6 < E < 6 and no solutions outside this ansKgjj p^nge. 



In many studies of the localization transitiona~Qll3~ll3 one 
looks for destruction of transmission at the band cen- 
tre, that is, only the behaviour of states at E = is 
considered. Since that is as far as possible from the 
unperturbed band edges, one is maximising the direct 
influence of quenched disorder which, coupled with the 



ever-present quantum interference effects, induces Ander- 
son localization (compare this with the purely quantum 
mechanism responsible for the very existence of a band 
of allowed states, and forbidden ones outside it, even in 
a homogeneous system). 

Further studies have investigated the E ^ region, 
attempting ta map the mobility edge trajectory in the 
E — W planelJuO. Along this latter line, reentrant be- 
haviour was found for certain distributions of site self- 
energies, in a narmw-jjange of energies above the unper- 
turbed band edgeLrlSo: at fixed E, as W decreases from 
some large initial value the system first displays an upper 
transition from insulator to metal at W — ► W£, much 
the same as at E = ; then, upon further reduction 
of disorder a second, lower transition takes place, back 
towards an insulating phase as W — > . The mech- 
anisms underlying each of these transitions argpqualita- 
tively understood, and fundamentally distinctQtl: while 
the onset of extended states as W — > W£ is due to dimin- 
ishing quantum interference, at W — » Wj~ it is the in- 
creasing overlap of (and consequent tunnelling between) 
localized states that eventually produces delocalized be- 
haviour. Reentrant regions were found for rectangular 
(P(£i) = constant, |ej| < W/2) and, to a lesser extent, 
for Gaussian (P(£i) oc exp(— ef/2W 2 )) distributions, and 
not at all for a Lorentzian one (P(£i) oc W/ [e| + IF 2 ]). 
This is consistent with the idea that, in order for local- 
ized states to coalesce into one extended wavefunction 
(the mechanism underlying the lower transition), their 
respective energies must not differ by much, thus en- 
abling their mutual resonance. The narrower the dis- 
tribution of self-energies (preferably with a sharp cutoff, 
as in the rectangular case), such situation is more likely 
to be found. For a Lorentzian shape, where already the 
second moment diverges, chances of overlap are negligi- 
ble. 

Though significant advances have beepjiiade recently 
regarding properties at the band centre£3il3, such as uni- 
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versality (with respect to varying disorder distributions) , 
corrections to scaling and accurate numerical estimates 
of critical quantities, a similar quantitative understand- 
ing of behaviour elsewhere on the phase diagram appears 
to be lacking. Here we set the model parameters in or- 
der to enable a study of the reentrant region, and probe 
both upper and lower transitions; our main goal is to find 
what universality class(es) they belong to. Our numeri- 
cal procedures and some concurring technical aspects are 
discussed in Section ||; the upper (high disorder) and 
lower (reentrant) transitions are analyzed respectively in 
Sections 
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given in Section 



conclusions and final remarks are 



II. MODEL DETAILS AND CALCULATIONAL 
METHOD 

We consider the Hamiltonian given in Eq. ([!]), on a 
simple cubic lattice, with a rectangular distribution of 
site self-energies: 



P{ei) = { 



constant - W/2 < e { < +W/2 
otherwise 



(2) 
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FIG. 1. Schematic phase diagram for rectangular disorder 
distribution (after Ref. 8, see their figure 1). Dashed line 
at E — 7 shows region investigated here, for which the lower 
and upper transitions are located at Wi — 8, W u — 15, re- 
spectively. 



For such distribution the mobility edge trajectory, as 
found in Ref. ||, is semi-quantitatively depicted in Fig- 
ure |l|. From inspection of the Figure, we decided for 
fixing the energy at E = 7 and scanning W. This way 
two well-separated transitions (respectively at W u ~ 8, 



Wi ~ 15) are expected to be present, so presumably one 
can focus on their individual properties one at a time. 

We have considered quasi- one-dimensional systems, 
"bars" of cross-section L x L sites, with periodic bound- 
ary conditions across, and very long length N. We have 
used L = 4, 6, 8, 10, 12 and 14; as explained below, the 
values of N were determined by numerical convergence 
criteria, 10 5 being typical. 

Following well-kpown procedures introduced by 
Pichard and Sarmafl (see also Refs. |6]JlO| for detailed de- 
scriptions), we considered the transfer matrix connecting 
matrix elements of the Hamiltonian, Eq. (|l|), on consec- 
utive cross-sections, and iterated it along the "infinite" 
direction, thus extracting the corresponding Lyapunov 
characteristic exponents. Of these, the most relevant 
here is the one with smallest modulus, which gives the 
largest localization length Xl(E, W) (slowest decay along 
the bar of transverse dimension L, at given W, for an 
electronic wave function with energy E). .— . 

According to finite-size scaling (FSS)B, the non- 
diverging quasi- one-dimensional lengths scale lin- 
early with L (apart from corrections to scaling, to 
be discussed below) at the critical point of the three- 
dimensional system. Therefore, plots of the scaled local- 
ization lengths Al = Xl/L (e.g. against W, for fixed 
E) for different values of L must cross, and they should 
do so twice if there are two distinct second-order tran- 
sitions. Leaving aside, for the moment, the question of 
whether this is consistent with the single-parameter the- 
ory of localization (see Ref. ||), one is provided with a 
rough-and-ready criterion to locate the expected critical 
points. 

In order to produce reliable numerical estimates from 
the conceptual framework described thus far, two points 
must be carefully considered. First, though the conver- 
gence of finite- N estimates of Lyapunov exponents to- 
wards— their asymptotic values is guaranteed in princi- 
pleMllZI, it may be hard to decide when it has taken place 
to a desired accuracy. As noted in Ref. ^, "for such very 
long strips the variation of A [as a function of N] from 
step to step is very small due to averaging, and this can 
lead to a false impression of convergence" . Accordingly, 
we have kept track of the variance among estimatesS, in 
the following way. Each 10 3 steps an estimate of A (tak- 
ing into account all iterations from the origin up to the 
sampling point) would be produced; first, twenty consec- 
utive estimates would be accumulated, and the average 
and RMS deviation among them calculated. At the end 
of each subsequent 10 3 steps, a new average and deviation 
would be calculated, of the most recent twenty estimates. 
Therefore, in this moving average procedure, one keeps 
track of fluctuations along the latest 20, 000 iterations. 
It was found that, in order to reach 0.1% fractional de- 
viation one typically needs 80, 000 — 130, 000 iterations 
(except for low disorder, W 6.5, where similar figures 
are achieved within ~ 50,000 steps or fewer). Strictly 
speaking, this is no guaranteee that the corresponding 
final estimate will be accurate within 0.1%, as a very 
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slow drift may still be masked underneath apparently 
random fluctuations. In Ref. [ll], 10 6 — 10 7 iterations at 
the band centre are quoted as necessary in order to give 
0.1% — 0.05% accuracy, although the latter term is not 
assigned a specific operational definition; in the closely 
related Green's function calculations of Ref. ||, 1% accu- 
racy is claimed for ~ 40, 000 — 70, 000 steps. Having this 
in mind, it seems that our current definition of 0.1% rel- 
ative fluctuation warrants perhaps 0.3% — 0.4% accuracy 
for the localization length, at best. We used 0.1% rela- 
tive fluctuations for all L = 4 — 10 data, and for selected 
L = 12 — 14 points; in the latter case, intermediate- 
ly data were collected with 0.5% — 1% fluctuations and 
checked for smooth variation against the more accurate 
ones. As shown below, this will not be the most impor- 
tant source of errors in the present case. 

Secondly, in general there may be systematic dis- 
tortions to the linear, single-scaling parameter picture, 
caused biLjionlinearity of scaling fields and/or irrelevant 
variables^. A detailed study of such effects for the An- 
derson transition at the band centre has been undertaken 
recently^. The result was that such corrections are of 
sizable importance in the case: estimates of critical quan- 
tities produced when they are taken into account tend to 
be more consistently reliable and accurate than if they 
are disregarded. Following the notation of Ref. O, near 
the critical disorder W c , with w = (W c — W)/W c , the 
scaled localization length A^ must vary as: 
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where \ is the relevant variable, ip is the leading irrel- 
evant variable, the exponent v characterizes the diver- 
gence of the thermodinamic (L — > oo) localization length 
via Ap q m w~ v , and y < is the leading irrelevant ex- 
ponentllaila. Since there is no transition for finite L, one 
assumes 
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further, since (for finite L) x and ip also are smooth func- 
tions of w, one expands both in power series as well; lin- 
earity is recovered if x( w ) — w ( as x(0) = must hold), 
ip(w) — ip (constant), and only the lowest powers are 
kept (nj = 1) in the Taylor series for A^: 



A, 



(5) 



Nonlinearity (both in relevant and irrelevant variables) 
and higher-order corrections associated to the irrelevant 
variable may be incorporated by keeping higher-order 
terms in the respective expansions. This is done at the 
expense of having a larger number of parameters to fit, 
and must be accompanied by the corresponding quality of 
numerical data. Here, we found that the quality of our 
data was generally consistent with the linearized form, 
Eq. (§). 



III. UPPER TRANSITION 

In Figure | we display the raw data for scaled local- 
ization lengths near the upper transition, whose approxi- 
mate location is estimated from Ref. |^ (see also Figure |l| 
above). Though data for each L behave somewhat ir- 
regularly, there is an unmistakeable crossing of curves, 
roughly at the centre of the diagram. 
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FIG. 2. Raw data for scaled localization lengths near the 
upper transition. Relative fluctuations, as defined in the text, 
are of order of symbol sizes or smaller. 

Following along the lines proposed in Ref. O, we at- 
tempted to fit assorted corrections to scaling. We found 
that, to be consistent with the range of variation and 
accuracy level of our data, the most sensible choice was 
(1) not to take nonlinearities of scaling variables into ac- 
count, and (2) keep only up to quadratic terms in the 
Taylor expansion of FifioL 1 '") in Eq. (||) above. Thus 
we plotted 
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corrected 
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ipL v (1 + ax + bx 2 ) 
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against x = (W — W u ) L 1 ^ . This gave us six parame- 
ters to fit: W u , v, ip, y, a and b. We started from reason- 
able guesses for some of these quantities (e.g. W u ~ 15, 
< 1); also, since this transition is driven by the same 
mechanism as that at the band centre, they might plau- 
sibly be in the same universality class, so v ~ 1.4 — 1.7 
seemed worth investigating (recall v — 1.57(2) from 
Ref. [ll]). With a few similar aditional assumptions, we 
checked the vaijiation of x 2 against changes in the fit- 



ting parameteraij; the plot displayed in Figure |3j is rep- 
resentative of the best sort thus obtained. Though its 
(unweighted) x 2 = 0.018 for 110 data does not mean 
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much by itself, it will be useful for comparison with the 
corresponding one at the lower transition, in Section [Tv| . 
Corrections to scaling are of order up to 5% for L = 4 and 
decrease steeply, owing to the large value of \y\ (similar 
to the picture at the band centre, see Table II of Rcf. 11 ), 
to around 0.2% for L = 14. The latter figure is close to 
the intrinsic numerical accuracy of raw data, so it is fair 
to say that corrections only play a significant role for the 
smaller sizes, say up to L = 8 — 10. Our plots turned out 
to be quite sensitive against changes in v of order 3 — 4% 
around the best-fitting value ~ 1.6; therefore, the initial 
hypothesis of universality, as incorporated in the fitting 
procedure, appears to be rather consistent. 
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FIG. 3. Scaling plot near the upper transition. Here, 
W u = 15.0, v = 1.61, i> = -0.21, y = -2.8, a = -0.21, 
b = -0.2 (see Eq. (g)). Lines are guides to the eye, starting 
at Ao = 0.568 on vertical axis. Key to symbols is the same as 
in Fig. |. 



On the other hand, the limiting amplitude, Ao = 
lhxiL->oo,W=W c given on scaling plots such as Figure || 
by the intersection with the vertical axis, was found to be 
remarkably robust: we quote Ao = 0.57(1) for all reason- 
ably behaved choices of parameters. This is in very good 
agreement with A = 0.576(2) , a value found to be uni- 
versal with respect-to the form of disorder distribution 
at the band centreEJ. 

TABLE I. Estimates of critical quantities at upper transi- 
tion. Data for the band centre from Ref. 11. Comparison of 
values for W c is not pertinent. 

W c A v y 

Present work 15.0(1) 0.57(1) 1.60(5) -2.8(5) 
Band centre 16.54(1) 0.576(2) 1.57(2) -2.8(5) 



Note also that Ao does not directly depend on the ini- 
tial choice of fitting parameters; instead, it is obtained 
at the end of the procedure, so the likelihood of this be- 
ing a biased result is smaller than e.g. in the case of 
v. Universal critical amplitudes play an important_ci 
in the statistical mechanics of critical phenomenaH 
so the present numerical agreement must be regarded as 
strong evidence in favour of universality of the transition 
(at least along the upper, non- reentrant, portion of the 
mobility edge, including the band centre). 

To our knowledge, the only existing calculation of ex- 
ponents along the mobility edge is that of Ref. g for 
Gaussian disorder (see their Table IV); however, their 
estimates are far too scattered, and no clear conclusion 
can be drawn from them. 

In Table Q we summarize our estimates for critical 
quantities at the upper transition; for comparison, cor- 
responding values at the band centre are quoted. Error 
bars for the present work's results are somewhat arbi- 
trary, but certainly encompass all well-behaved scaling 
plots. The exact coincidence of the estimate for y with 
that of Ref. [ll] must be regarded as fortuitous. The over- 
all conclusion is that both transitions are in the same 
universality class. 



IV. LOWER (REENTRANT) TRANSITION 

In Figure |] are displayed the raw data for scaled local- 
ization lengths near the lower transition, whose approxi- 
mate location is again estimated from Ref. g. 
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FIG. 4. Raw data for scaled localization lengths near the 
lower transition. Key to symbols is the same as in Fig. ^. 
Inset: raw data for L = 6 (full line) and L — 12 (dashed line), 
2 < W < 16.5. 
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Data here behave rather differently from the upper 
transition. Non-monotonicity has been found earlier in 
this region; see data of Ref. ^ for E — 6.6, 6.8, W = 8, 
L = 6 — 10. A proper perspective is gained by looking at a 
broader range of IF-variation, including both transitions: 
see the inset of Figure ||, where only L — 6 and 12 data 
are displayed, for clarity. The maxima shown in the main 
diagram are in fact shoulders, attributable to irrelevant 
scaling variables, much in the same way as the corrections 
discussed in Ref. |ll| and Section III above. They must 
vanish for increasing L, and are distinct from the bona 
fide maxima, expected to be exhibited by finite-size data 
deep inside the conducting phase, somewhere between 
the two transitions. Indeed, the inset clearly shows that 
for L = 6 the former effect is stronger than the latter, 
while for L = 12 the trend has very much been reversed. 
For L = 4 (not shown) shoulder and maximum (if any) 
merge, owing to finite-size broadening of the former. 

The magnitude of distortions thus displayed makes 
it hard to account for them within a Taylor expansion 
framework, with a realistic number of parameters. At- 
tempts to subtract the shoulders by postulating an ad 
hoc gaussian shape were not successful. We then decided 
to discard L = 4 — 8 data, and try a fitting procedure 



similar to that given in Eq. 



j) above, using only data for 



L = 10, 12, and 14. Figures shows the results of the best 
fit, with a set of parameters for which the (unweighted) 
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FIG. 5. Scaling plot near the lower transition. Here, 
Wi = 7.9, v = 1.65, tp = 4.5, y = -3.2, a = -2.07, b = (see 
Eq. ([]))• Lines are guides to the eye, starting at Ao = 0.54 
on vertical axis. Key to symbols is the same as in Fig. ^; only 
data for L = 10 — 14 are used. 



One sees that the shoulder of L = 10 data has been re- 



duced, but not totally eliminated; also, the overall quality 
of data collapse is inferior to, and the amplitudes of cor- 
rections larger than, the corresponding ones at the upper 
transition. However, the values v — 1.65 and Ao = 0.54 
used in Fig. || are within 6% of those at the band centre. 
Bearing in mind the numerous assumptions made along 
the way, assigning error bars to these numbers would 
be somewhat risky. We have found that plots using v 
just about anywhere in the 1.4 — 1.7 range do not seem 
obviously much worse than the one displayed above; as 
regards A , anything between 0.5 and 0.6 seems possi- 
ble, provided that one keeps Wi in the range 7.8 — 8.2. 
Therefore, the proper conclusion is that our data are not 
inconsistent with the lower transition being in the same 
universality class as that at the band centre (and, pre- 
sumably, along the upper portion of the mobility edge). 



V. CONCLUSIONS 

We have found strong evidence that the transition 
along the upper portion of the mobility edge is in the 
same universality class as that at the band centre. We 
quote v = 1.60(5) and A = 0.57(1), to be compared re- 
spectively with v = 1.57(2) and A = 0.576(2) of Ref. 0. 

Our results for the lower (reentrant) transition are 
more mixed, though they certainly are compatible with 
universality: we give v — 1.55(15) and Ao = 0.55(5), 
from a qualitative visual analysis of trial scaling plots. 
One likely source for this would be the significant pres- 
ence of finite-size corrections, linked to irrelevant scaling 
variables and/or nonlinearities in the scaling fields; so 
far, we have not been able to account for the bulk of 
such effects within a Taylor-expansion picture similar to 
that successfully employed, both at the band centre in 
Ref. [ll] and at the upper transition in the present work. 

Another, more fundamental, reason for such irregulari- 
ties would be the failure of one-parameter scaling theory 
to describe the lower transition at all. It is argued in 
Ref. H that "the critical behaviour should be independent 
of whether the energy E, the disorder W, or some com- 
bination of both, is taken as the scaling variable ... this 
is not generally true [e.g., in the present case] because 
the system shows two metal-nonmetal transitions when 
the disorder is changed for a fixed energy" . It is also true 
that the physical mechanisms driving the lower transi- 
tion (quantum tunnelling between tail states) and the 
upper one (quantum interference between essentially de- 
localized wave functions) are, in principle, distinct. Thus 
perhaps the one relevant parameter that represents the 
latter physical effect cannot do so for the former. 

On the other hand, reasonably good scaling plots have 
been found here, by implicitly assuming the validity of 
one-parameter scaling, with corresponding numerical val- 
ues of critical quantities which are not inconsistent with 
universality. This hints that the picture along both the 
upper and reentrant parts of the mobility edge may, in 
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fact, be rather simple: critical behaviour would be dom- 
inated by a single fixed point, located at E = 0. 

More work, numerical and theoretical, is clearly nec- 
essary in order to settle the issue. It is also expected 
that experimentalists could be attracted into the discus- 
sion; metallic samples with suitably-located Fermi ener- 
gies might conceivably be made to scan reentrant regions, 
by varying dopant concentrations. 
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